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Abstract 

We study material lines and passive vectors in a model of turbulent flow at infinite- 
Reynolds number, the Kraichnan-Kazantsev ensemble of velocities that are white-noise in 
time and rough (Holder continuous) in space. It is argued that the phenomenon of "spon- 
taneous stochasticity" generalizes to material lines and that conservation of circulations 
generalizes to a "martingale property" of the stochastic process of lines. 
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The evolution of material lines and surfaces passively carried by turbulent flow has long 
been a subject of interest [1]. This is motivated in part by questions surrounding dynamically 
relevant objects, such as vortex lines O O S] and magnetic field-lines [SJ [U], which have been 
argued to behave similar to material lines at high Reynolds numbers. However, in that limit, the 
turbulent velocity field is no longer differentiable in space but only Holder continuous [3 El E] . 
Observations from experiments and simulations suggest that material objects advected by such a 
rough velocity become fractal, with a Hausdorff dimension strictly greater than their topological 
dimension [TT] [T^ [THl IT^ IT^ . This poses a difficulty to the view that vortex lines behave 
as material lines — a consequence of the Kelvin-Helmholtz theorem [T7] — since circulations 
are a priori not defined for non-rectifiable loops. It has recently been argued that the Kelvin 
theorem in fact breaks down in turbulent flows, in the sense that the circulation is not strictly 
conserved for every loop |18^ I19j . Similar breakdown of Alfven's theorem on magnetic-flux 
conservation [20j is expected in plasma turbulence at high magnetic Reynolds numbers |21j . 

These questions have been sharpened by recent work on the Kraichnan model of advection 
by a Gaussian random velocity field that is delta-correlated in time [22|. A novel phenomenon 
has been discovered there called spontaneous stochasticity: Lagrangian particle trajectories for 
a non-Lipschitz advecting velocity are non-unique and split to form a random process in path- 
space for a fixed velocity realization [231 1211 [251 ISSl IIZI [281 [29]. This phenomenon raises many 
fundamental questions, including whether material objects such as lines and surfaces can even 
exist in the limit of infinite Reynolds number. It is the purpose of this Letter to outline a 
new approach to the evolution of such geometric objects in the Kraichnan model. We focus on 
material lines and passive vectors, which are dual objects in the same sense as material particles 
and passive scalar fields [30]. In particular, we shall sketch the proof of a "martingale property" 
that has previously been proposed [18] as a generalization of the conservation of circulations 
for a rough velocity field. 

We consider stochastic flows |31j on a d-dimensional manifold A4 driven by Brownian vector 
fields that are not Lipschitz regular in space. To simplify the presentation, we use Euclidean 
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space A4 = M.'^ OT the torus = T*^ to illustrate the main ideas. More precisely, u(x, t) is a 
Gaussian random vector field, with mean u(x, t) and fluctuating part u(x, t) with covariance 

(Si(x,t)Sj(x',t')) = Ai(x,x';t)5(t-0- (1) 

for X, x' € A^. We are mainly interested in the case that u(x, = and u(x, t) = u(x, t) is a 
homogeneous random field, with D(x,x';i) = D(x — x',i). The quantity 

A(x, x'; t) = tr[D(x, x; t) + D(x', x'; t) - 2D(x, x'; t)] (2) 

is (||u(x) — u(x')lp), the mean of the Euclidean norm squared, for a random field u(x) with 
covariance D(x, x';t). The case of greatest interest to us has A(x, x';t) oc ||x — x'p" with 
< a < 1. In that case, u(x) is Holder continuous with exponent a at every point in space. 

We consider oriented lines (1-cells) given parametrically as continuous, one-to-one maps 
C : [0, 1] ^ A^. A material line satisfies 

{d/dt)C{a, t) = u(C((j, t),ot). (3) 

for a G [0, 1] and t € M. The circle "o" means that we interpret equation ([3]) in the Stratonovich 
sense. The (forward) Ito equation (d/dt)C{(T,t) = u(C{a,t),t) equivalent to ([3]) has the mean 
changed to u*(x,f) = Ui(x, t) + (l/2)(5/5x'=)Afc(x, x'; f)|x'=x ([31], section 3.4.) If X = or 
T'^ and if u(x,f) is a homogeneous random field, then the Ito and Stratonovich interpretations of 
equation ([3]) are equivalent. Now let Pu[C')i] denote the conditional probability distribution of 
lines for a fixed velocity realization u. This distribution satisfies a stochastic Liouville equation: 

{d/dt)P^[C,t] = - [\a—^{u,{C{a),ot)P^[C,t]). (4) 
Jo dCi{a) 

Equation is a direct consequence of equation ([3]) and must also be interpreted in the 
Stratonovich sense. It is formally equivalent to the Ito equation: 

{d/dt)P^[C,t] = - [\a—^{[u*{C{a),t)+Ui{C{a),t)]Pu[C,t]) 
Jo oCi[a) 

^\[/'jy'K^£(W) (A,(C(,),a(.');t)Pu[C,tl) . (5) 
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Averaging equation ([5]) over the Gaussian ensemble of velocities u yields a functional Fokker- 
Planck equation for distributions in the space of free lines C on the manifold Ai : 

(dldt)P\C,t\ = - j\.^A- ([ul(C{a).t)P[C.t]) 

^U'-l'-^'mScW) (^«(C(<')-C'(-');t)P[c.*l) . (0) 

The first term on the righthand side represents a drift with the mean velocity u* and the 
second term represents a diffusion arising from the velocity covariance D. Similar diffusions 
on the path- and loop-spaces of a manifold A4 have been much studied, motivated in part by 
questions from quantum field theory |32[ [55] . 

The above considerations are rigorously justifiable for the case of a Lipschitz velocity with 
a = 1 but are only formal when a < 1. A more careful (and also more physically realistic) 
approach in the latter case is to replace the advecting velocity u with a "coarse-grained" or 
smoothed velocity ux = (px* u, by convolution with a smooth filter kernel 93A(r) = X^'^ip{r / X). 
The length-scale A can be interpreted as a mathematical representation of the viscous cutoff in 
a true turbulent velocity field . The exact solution of the Liouville equation (jH) for such 

a smoothed velocity is 

PuMC,t\Co,to) = S{C - $'^'{Co))dC, (7) 

1 1' 

with initial condition Pu^{dC,tQ\CQ,tQ) = 5{C — Co)dC. Here : ^ is the stochastic 
flow of diffeomorphisms generated by the smoothed velocity- field ux (131], section 4.6). Despite 
([7]), a nontrivial diffusion process in line-space can be obtained if the limit A — > is taken 
appropriately. Consider a "nice" distribution Gp{dC) which is supported on lines entirely 
contained in the ball B{0,p) of radius p at the origin and take the weak limit 

hm hm [ GpidC'o) [ Pu,{dC,t\Co + C^,to)^{C) = [ Pu{dC,t\Co,to)^{C) (8) 

for bounded, continuous functionals ^'(C) and t > tQ. In the "weakly compressible regime" 
[24:\ — and, in particular, for a divergence-free velocity field satisfying V-u = — this 
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limit should yield a non-degenerate diffusion. This is a generalization of the phenomenon 
of spontaneous stochasticity to the turbulent advection of lines, with initial line Cq at time to 
splitting into a random ensemble of lines C at time t. 

As for the case of smooth advection, an unconditional diffusion satisfying equation ^ 
may be obtained by averaging over the velocity u. The instantaneous realizations C of this 
diffusion process should be fractal objects when the advecting velocity is Holder continuous 
with exponent < a < 1 and rigorous estimates of their Hausdorff dimensions would be of 
much interest. These questions may also be addressed numerically using Lagrangian Monte 
Carlo techniques \34: \ 135 1 [36]. such a study, the material line C{t) would be represented by a 
discrete approximation CN{t) constructed from A^-|- 1 Lagrangian particles Xa(t), a = 0, : 

CN{cT,t) = (1 - 6lAr(cr))Xaj^(^)(t) + 6lAr(cr)Xaj^(^)+l(t). (11) 

Here aN{a) = [Na] with [x] the greatest integer less than or equal to x (modulo N for loops) 

and 9n{(^) = {^(^) where {x) = x — [x] is the fractional part of x. Thus, (fTT|) corresponds to a 

piecewise-linear curve with linear segments connecting the successive Lagrangian particles. So 

long as 5N{t) = max^ |xa(t)— Xa_|_i(t)| < A, the discrete approximation CAr(t) represents well the 

material line C{t) and the approximation becomes better as ^ oo and 5N{t) <C A. However, 

the same is not true in the opposite limit, where A <C ^^{t). The phenomenon of "spontaneous 

stochasticity" for a rough velocity field makes it very doubtful that material lines even exist if 

the limit A — > is taken before evolving in time and an initial line then presumably "explodes" 
^We note that the same diffusion process should also be obtained by using Duhamel's formula to solve equation 
((5]) as the Ito integral 

Pu[C, t] = S*o{t, to)P[C, to] - f dt' S*{t, t') [ da—^ {[ui{C{a),t')]P4C, t']) , (9) 

Jto Jo SCi{a) 

where So{t,t') = Texp [j^*^ dt'A*{t')^ and 

is the Fokker-Planck operator of the (mean) diffusion in line-space. Equation ([9]) can be solved iteratively to 
generate a representation Pu[C, t] = S^{t,to)P[C,to] as a Wiener chaos expansion in white-noise u; cf. [281 129] . 
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into a disconnected cloud of particles at any time t > 0. Thus, the velocity smoothing in ([8|) 
appears to be necessary to define appropriately a line-diffusion for a rough (Holder) velocity. 
Alternatively, a stochastic regularization might be employed that adds a white-noise KdW{t) 
to the evolution equation of Lagrangian particles |30j . 

We now turn to the dual problem of a passive vector (1-form) A advected by the velocity 
u = u -|- u: 

(9tA(x, t) + (u(x, ot). V)A(x, t) + (Vu(x, ot))A(x, t) = 0, (12) 

This stochastic equation is interpreted again in the Stratonovich sense. Equation (|12|) for d = 3 
is equivalent by vector calculus identities to dtA + V(u-A) — ux(VxA) = 0. The latter has 
the form of Ohm's law, 

E + uxB = ??J, (13) 

in the ideal limit of zero resistivity (r/ = 0) for an electric field E = —dtA — V(u'A) and 
magnetic field B = Vx A given by a vector potential A, with the electric current J = VxB. 
Taking the curl of (jlSp yields an induction equation S^B = Vx(uxB)-|-ryAB for the magnetic 
field. With this interpretation, the passive vector equation was introduced by Kazantsev |37J 
as a soluble model of the kinematic dynamo. (See also \38\ [39} HOj.) The "circulation" (or 
"holonomy") of A along C is defined in Lagrangian form as 

^LiC,t) = [ A(t).dx, (14) 

Jc{t) 

where C{t) is the material line advected by u(t) which started as line C at the initial time 
t = to- Conservation of "circulation", {d/dt)^L{C,t) = 0, follows formally from (jl2p for any 
space dimension d > 1. It is rigorously true for the case of a smooth advecting velocity with 
a = 1 ([31], section 4.9). If C is a closed loop (1-cycle), then the line-integral (fH|) represents 
gauge-invariant magnetic flux and the conservation law corresponds to Alfven's theorem |20| . 
We now consider the generalization of this result for a < 1. For this purpose, it is useful to 
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reformulate the passive vector equation ()12p as a passive scalar in line-space: 



dt^E{C,t)+ [ da [ui{C{a),t)+u,{C{a),ot)] 
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^E{C,t) = 0. 



(15) 



Jo 



In this equation, $£;(C, i) = A{t)-dx is the Eulerian circulation of A along a fixed (non- 



advected) line C. An exactly analogous reformulation of the incompressible Euler equation 
(as an active scalar in loop-space) was advanced some time ago by Migdal ^41j. Note that 
conservation of circulations is just the formal solution of (jlSp by the method of characteristics. 
We shall take the equation (|15p as our primitive formulation of the passive vector; one of 
the immediate advantages is that we can avoid (for the moment) the question how to define 
line-integrals over fractal lines. We then convert (jlSp to Ito formulation: 



This stochastic equation is solved by the method of LeJan and Raimond [28\ [29] (cf. footnote 
7^1), writing it as a (backward) Ito integral and iterating to obtain ^EiC, t) = Su{t, t')^E{C, t'), 
where the Markov operator semi- group Su{t,t') is defined by a Wiener chaos expansion. More 
intuitively, this solution is expressed as 



in terms of the turbulent diffusion of lines (backward in time). We see that the circulations 
are not conserved, except on average. This is precisely the "martingale property" that was 
conjectured (for solutions of incompressible Euler equations) in [18]. Note that this property 
imposes an irreversible arrow of time, since Eulerian circulations are given as averages over 
their past values, not future ones. This "generalized Alfven theorem" should be related to 
dynamo action in the Kazantsev model |37l I38[ \39 \ 140] . In the physical context of the dynamo, 
there is an additional resistive term rj §^dxjdij^^-j^<^E{C,t) on the righthand side of (fT5]) . 
where 5 / daij (x) is the "area derivative" in the loop calculus of Migdal [H] . This term breaks 
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(16) 




t' < t, 



(17) 



time-reversal symmetry and should select the backward-martingale solution (jl7p in the ideal 
limit r/ — > 0. Of course, it cannot be ruled out a priori that the r] ^ limit of the resistive 
regularization and the A ^ limit for regularized velocity, as in ([SD, shall yield distinct weak 
solutions of the loop-equation (|15|) . as occurs in the intermediate compressibility regime of the 
passive scalar problem [5H [251 123 EB] . 

These results can be generalized to turbulent diffusion processes of higher-dimensional ma- 
terial objects, /c-dimensional oriented submanifolds of A4 or k-cells C^{t). The dual object is 
the passive k-form lo^, which satisfies (in Stratonovich sense) 

dtuj'' + Luuj'' = (18) 

with Lu the Lie-derivative along the vector field u ([31], section 4.9). This equation is formally 
equivalent to conservation of the integral invariants 

I{C\t)=f u\t) (19) 
Jc^(t) 

for any fc-cell C^{t) comoving with u f42]. Then A; = is the passive scalar, k = 1 the passive 
vector and k = d the passive density |2l]. A theory similar to that developed here for k = 1 
applies for any integer k. A unified approach to all these results is to consider directly the 
turbulent diffusion of the Lagrangian flow maps , which satisfy the stochastic equation 

(dM)^*'*'(a) = u(^*'*'(a),ot), ^*''*'(a) = a. (20) 

In this framework one can derive for the distribution Pu[^,^] on maps exact analogues of the 
Liouville equation, in Stratonovich form ^ or Ito form ([5|). It is natural to formulate the 
problem as an infinite-dimensional diffusion in the Hilbert space n = L2(X,M^) . It is known 
for the cases M = R'^ or T*^ that the semigroup S{M) of Borel volume-preserving maps is a 
closed subset of this Hilbert space, and that the group G{A4) of volume-preserving diffeomor- 
phisms is dense in S{A4) for the L^-topology [43j. This construction is a close analogue of the 
"generalized Euler fiows" of Brenier, but for the Cauchy initial-value problem. 
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To summarize: We have outlined an approach to the study of material lines in a model 
of turbulent flow at infinite Reynolds number and to the dual problem of a passive vector in 
the same flow. The main conclusions are (1) that a non-degenerate diffusion should exist for 
material lines, generalizing the phenomenon of "spontaneous stochasticity" of material points, 
and (2) that the Kelvin/Alfven theorem on conservation of circulations should generalize to a 
"martingale property". Although the approach sketched here depends heavily on the white- 
noise character of the velocity field in time, we expect that similar results hold for more realistic 
velocity ensembles with the crucial property that realizations are rough (Holder continuous) in 
space. See |181I21| for related rigorous results on the solutions of incompressible fluid equations. 
The two properties discussed in the context of this model problem should be an essential feature 
of real fluid turbulence in the high Reynolds number limit. 
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